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The Kibble-Zurek mechanism (KZM) captures the key physics in the non-equilibrium dynamics of 
second-order phase transitions, and accurately predict the density of the topological defects formed 
in this process. However, despite much effort, the veracity of the central prediction of KZM, i.e., 
the scaling of the density production and the transit rate, is still an open question. Here, we 
performed an experiment, based on a nine-stage optical interferometer with an overall fidelity up 
to 0.975±0.008, that directly supports the central prediction of KZM in quantum non-equilibrium 
dynamics. In addition, our work has significantly upgraded the number of stages of the optical 
interferometer to nine with a high fidelity, this technique can also help to push forward the linear 
optical quantum simulation and computation. 



In the early universe after the " Big Bang" , cosmolog- 
ical phase transitions occurred with the expansion and 
cooling of the universe, and the symmetry of the vac- 
uum was broken. The new vacuums were chosen lo- 
cally, within space-like regions, resulting in topological 
defects [J [5]. The initial density of the topological defect 
is extremely interesting, and a rough limit of this density 
can be estimated by the light-cone causality; however, 
the exact density is not easy to determine. Zurek sug- 
gested that this cosmological mechanism can be observed 
in condensed matter systems in a laboratory [3™S]. For 
example, a pressure quench drives liquid 4 He from a nor- 
mal phase to a superfluid phase at a finite rate, which 
leaves behind vortex lines. 

In a condensed matter system, the speed limit of the 
light is less useful to estimate the density of defects. How- 
ever, the density of topological defects can be predicted 
for second order phase transitions due to the divergence 
of the relaxation time r (which characterizes the time 
required for the order parameter to relax to its equilib- 
rium value when the parameter has been perturbed) and 
the healing length (which characterizes the length over 
which the order parameter will return to the equilibrium 
value when disturbed) near the critical point. As a re- 
sult of the divergence, every such transition, traversed at 
a finite rate, is inevitably a non-equilibrium dynamical 
process. The whole system can not catch up, and the 
symmetry will be broken with some topological defects 
[3H5]. Therefore, the Kibble-Zurek mechanism (KZM) 
provides a theoretical framework with which to describe 
the non-equilibrium dynamics of the symmetry broken in 
the second order transition [THl]. 

The central prediction of the KZM is that the den- 
sity of the topological defects in the second order phase 
transition should scale with the transit rate [3H5]. To 
be more specific, consider the phase transition in liquid 
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4 He driven by the pressure (denoted here by t) and sup- 
pose that the critical point is located at t = 0. The 
entire dynamic process can be divided into three parts, 
the adiabatic, impulse and adiabatic regions(which are 
shown in Figjlja)), according the relaxation time scale 
r. Between —t and t, the system can not adjust in time 
to follow the new thermodynamic conditions and is ef- 
fectively frozen. The freeze-out time t is determined by 
the Zurek equation [3H5]: r(t) = t. The relaxation time 
r with a general parameter t can be determined experi- 
mentally as t = To/\e\ for the case of liquid 4 He (where 
To is a constant, and e with e(t = 0) = is called the 
relative temperature, which is used to characterize the 
distance between the system and the critical point). In 
addition, e can be related to the transit rate tq in the 
liquid 4 He case as e(t) — t/rq. Therefore, the freeze- 
out time can be dcterminged with the Zurek equation as 
t = y/TQTQ. As analyzed in [3H5] , the density of the topo- 
logical defects produced in the non-equilibrium dynamics 
of the transition can be predicted with the knowledge of 
t without solving the dynamical equations. 

For the ubiquity of the second order transition and 
non-equilibrium dynamics, it is very interesting to verify 
the KZM in the laboratory. Though many efforts have 
been made[oT-fl~5]. the central prediction of the KZM has 
still not been observed clearly [I5J. Because the KZM is 
a framework for non-equilibrium dynamics and the an- 
ticrossing between the ground state and the first excita- 
tion state of a generic second order phase transition[T7], 
Damski [TS] presents a simple scheme to support the 
KZM in a quantum non-equilibrium dynamical system 
- the Landau-Zcncr (LZ) model with the time dependent 
Hamiltonian H = (uj a x + At<r z ) /2 [T8H2U] . By introduc- 
ing correspondence between the LZ model and the second 
order phase transition, such as, tq := ujq/A, tq := 1/loq 
and the density of defects D n := |(\l/|l)| 2 , and based on 
the framework of the KZM, Damski provides a direct pre- 
diction of the transition probability long after the critical 
point, which corresponds to the density of the topological 
defects in the second order phase transition (The details 
are shown in FigjT] and the procedure for comparison is 
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FIG. 1. (color online) The comparison of the features in the KZM and the quantum non-equilibrium dynamics of the LZ model, 
(a) shows the relaxation time scale r in the KZM, and (b) shows the inverse of the energy gap in the LZ model, (c) gives the 
structure of the energy levels in the LZ model. 



given in the supplementary information). Suppose the 
dynamics start at t = 0, the transition probability at 
t > under the KZM should then be 

V n = 0.5(1 - y/l-2/V(x a )), (1) 

where x a — a^r and V(x a ) = x 2 a + x a ^/x^ + 4 + 2 
(see the supplementary information and [TS] ) . Eq[l] is 
a direct result of Zurek's equation and reflects the key 
quantitative aspects of KZM, i.e., the scaling rule, in this 
quantum dynamical process. 

To verify Eq[TJ which is predicted by the KZM, we 
need to realize a quantum dynamical process, which is 
driven by a time-dependent Hamiltonian. The dynami- 
cal process is not easily implemented in linear optics. The 
traditional methods reported as quantum simulations in 
linear optics [21] . which require the prior information of 
the state and can be viewed as a state preparation pro- 
cess, can not be used in this situation. Although the 
information of the state of the simple LZ model can be 
calculated by solving the dynamical equations, the tra- 
ditional method is not a dynamical process and can not 
be used to verify the KZM. In the dynamical process, 
we only have the information of the initial state and the 
time-dependent Hamiltonian (or the corresponding uni- 
tary operator), and have no prior information about the 
middle state and the final state. We first experimentally 
implement this dynamical process in linear optics. We 
divide the dynamical process into many segments due to 
the difficulty of tuning the Hamiltonian continuously; in 
each segment (a step in the experiment), the Hamiltonian 
has a fixed parameter according to a different time. This 
process is a truly quantum simulation in linear optics 
and our method has been used to realize an imaginary 
time evolution and demonstrate the algorithmic quan- 
tum cooling [22] • Here, we extend the process to a real 
time evolution with more stages. With this method, we 
observe the evolution of the LZ model, and by comparing 
the experimentally measured transition probability and 
the prediction of Eq[l] the first experimental observation 
of the KZM in a non-equilibrium quantum dynamical 
system is presented. 



The basic module in our experiment, shown in 
Figj2ja), is a polarization based interferometer, known 
as the Mach-Zender interferometer (MZI), with some lo- 
cal operations that are realized by wave plates inside and 
outside of the device. The MZI is one of the most ba- 
sic elements in a linear optical quantum computation. 
It has been noted that, in principle, the quantum logic 
in a quantum computation can be simulated universally 
by an MZI [23]. Figurej^b) shows the lo gic diagram of 
this structure, and the overall operation is then given 
by ^= (I — i<j x sin 4a cos 27 — ia z cos 4cv cos (27)) where 

0(7) represents the angle of the optical axis of the HWPs 
(QWPs) and a x and er z are two Pauli operators. In dig- 
ital quantum simulators, to obtain the final state of a 
quantum system governed by a time-dependent Hamilto- 
nian H sys (t) for a time interval T, we divide the entire 
process into n steps. In each step, with a time inter- 
val t = T/n, H sys (t) will not change significantly and 
can be replaced with a time-independent Hamiltonian 
H(tk), where k is the index of the step. The time tk can 
be at any time position of the step k. For convenience, 
it is chosen here at the middle of the step. For a ini- 
tial state \<j>), the final state can then be approximated 
by life e~ lhHtytk ^ T \(f>). The only remaining task is to real- 
ize each step's evolution and, step by step, to then obtain 
the approximate final state. For each step's evolution, we 
can replace g-^^wr ( approximately, with I — ihH (ifc)r , 
where / is the identity operator. For the Hamiltonian in 
the LZ model, the approximate local time evolution op- 
erator is Uk = I — i(T x 6 — io'zOek, where 9 = and 
£fe = ~ k -. Comparing this operator with the operator 
formula given above, we can obtain the relevant value of 
a (7) and adjust the Hamiltonian in each step by chang- 
ing the angles of the wave plates. 

Due to the intrinsic phase instability in the standard 
MZI, it is impossible to build an interferometer with com- 
plicated structures and multiple stages. That is to say, 
we can not extract any practical quantum logic based 
on the standard MZI. Scientists have developed an in- 
creasing number of schemes to overcome this disadvan- 
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FIG. 2. (color online) (a) A standard MZI, whose logic diagram is shown in (b). (c) provides the experimental details. FC 
stands for fiber collector. 



tage [HJ [35] . Thus far, these methods have provided no 
more than five stages. Here, shown in FigureMc), we up- 
grade the number of stages in optical interferometer to 
nine based on another alternative device, the Sagnac-type 
interferometer, which has been used to realize a partial 
measurement [26] and has been developed to three stages 
to demonstrate the algorithmic quantum cooling [22] . In 
our experiment, the visibility of every stage is over 1000:1 
and the overall fidelity is measured to be 0.975±0.008, 
which is the highest known value that provides sufficient 
intrinsic stability to perform small-scale quantum com- 
putations. The characterization of nine-stage interferom- 
eter is given in detail in the supplementary section. 

The main purpose of this work is to show the Kibble- 
Zurek mechanism in the non-equilibrium dynamics of the 
LZ model. First, we choose the parameters A = 5 and 
t = 2 and perform a simulation of the evolution of the 
LZ model (The results are shown in the supplementary 
information). Due to the limitation of the number of 
stages, we have clearly observed one and a half periods 
in the evolution of the LZ model. Then, to achieve the 
main goal, we change the parameter coq from 0.05 to 2.5, 
which enables us to observe the LZ dynamics from the 
near equilibrium case to the non-equilibrium case. We 
compared the density of the H polarization after the 
evolution by the nine-stage interferometer and the the- 



oretical prediction by Eq[T] these values are shown in 
Figure[3] When the value of tq/tq is small, which is 
the case in the non-equilibrium region, we find that the 
density of the upper lever measured in the experiment is 
closely equivalent to the density of the topological defects 
predicted by the KZM (shown by blue dash-dotted line) . 
The green line gives the numerical solution of our experi- 
mental system, which deviates from the KZM prediction 
in the large region of tq/tq. There are two reasons for 
this deviation. First, the limited number of stages in 
our experiment (which we believe to be the maximum 
number of stages achieved thus far) means the evolution 
has not completely finished (shown in the supplementary 
section); second, the system has deviated from the non- 
equilibrium region, which can be confirmed by the exact 
analytical solution of the LZ model, shown by the pink 
dashed line in Figure (3j 

In conclusion, a nine-stage interferometer has been 
built with a high visibility and intrinsic stability. By us- 
ing this device, we observe the dynamical evolution of LZ 
model. As the main objective of this work, we change the 
quench time scale in the LZ model and measure the final 
density of the upper level which corresponds to the den- 
sity of the topological defects in the KZM. We then pro- 
vide an experimental confirmation that the KZM emerges 
from the simplest quantum non-equilibrium dynamics. 
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FIG. 3. (Color online) The experimental results compared 
with the theoretical prediction. The blue dash-dotted line is 
the theory prediction from the Kibble-Zurek mechanism with 
a = 0.8. The experimental results are shown by the points, 
and the green line shows the numerical solution of our ex- 
perimental system. The pink dashed line gives the analytical 
solution with the same parameters. 



We hope that these results will be helpful in understand- 
ing the KZM in quantum non-equilibrium dynamical sys- 
tem and in revising the simplest system from another 
perspective. Our platform will be helpful in quantum 
walks and quantum computations. When we completed 
our experiment and prepared this manuscript, we noticed 
that K. Pyka et al, [27], reported an experimental result 
to demonstrate the scaling of the number of topological 
defects with the transition rate in ion trap. 



METHODS 

The experimental setup is shown in Figure{2|c). Co- 
herent pulses from a Ti:sapphire laser with the center 
wavelength mode locked at 800 nm (with a 130 fs pulse 
width and a 76 MHz repetition rate) are adequately at- 
tenuated to the single-photon level, which is less than 
0.007 photons in each pulse, and coupled to a nine-stage 
polarization based interferometer via a single mode fiber 
(SMF). This resource has been employed to demonstrate 
quantum walks [28] and has also been widely used in 
practical quantum key distributions |29j . The initial po- 
larization state is prepared by using a polarized beam 
splitter (PBS) followed by a half wave plate (HWP). Each 
step of the evolution process is realized by a module com- 
posed of a polarization-dependent Sagnac interferometer. 
By using two quarter-wave plates (QWPs) independently 
in each arm and two HWPs (one in front of and the other 
after the interferometer), we can adjust the evolution pa- 
rameters; this process is described in detail in the main 
text. A phase compensator (PC) is inserted in one arm 
to correct the phase error and there is an adjustable slit 
in the other arm for calibration. Aim long SMF is used 



as the coupler between two independent interferometers 
and also acts as a spatial filter. When the evolution is 
finished, the photons are collected by an SMF and guided 
to the polarization analysis module, which is consists of 
an HWP, a Wollaston prism and two single photon de- 
tectors (SPDs). One interference filter (IF) with a 1 nm 
band width and an 800 nm center is introduced before the 
Wollaston prism to increase the coherence length. Two 
QWPs, one after the state preparation and the other be- 
fore the polarization analysis module, are inserted to per- 
form the process tomography. 

About the errors: In our work, the errors are pri- 
marily due to three factors: The first error is the errors 
of the wave plates. There are more than 30 wave plates 
in the setup, and each one has an angle error of approx- 
imately 0.1 degree: the total error can be simulated nu- 
merically and is found to be 0.011 by the Monte Carlo 
method. The second error is the shot noise. In our ex- 
periment, the total number of photons is over 5 x 10 5 : 
according to the numerical calculation, the error due to 
the shot noise is much smaller than the angle error of the 
wave plates. The last error is due to the power fluctua- 
tion of the photon source: here, we employ two detectors 
to detect all of the outputs and overcome this disadvan- 
tage by renormalization. 
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SUPPLEMENT FOR QUANTITATIVE 
VERIFICATION OF THE KIBBLE-ZUREK 
MECHANISM IN QUANTUM 
NON-EQUILIBRIUM DYNAMICS 

A. Landau-Zener model as the simplest version 
supporting the Kibble-Zurek Mechanism 

The LZ model is a well understood quantum dynam- 
ical model for the relatively simple form of the time- 
dependent Hamiltonian Hlz = 0.5(Ata z -\-u)Q<y x ), which 
has an analytical solution in theory. In this model, the 
so called LZ transition acts as a result of the tuning, and 
the transition probability depends on both the energy 
gap in the critical point, ujq, and the velocity of the tun- 
ing of the LZ Hamiltonian, A. The former is decided by 
the quantum system itself; increasing the gap decreases 
the probability. The latter is related to the speed of the 
tuning of the external field. These two parameters deter- 
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mine the instant energy gap, a/^o + (At) 2 . According 
to the adiabatic theorem, when the inverse of the gap is 
small enough, the system can catch up with the changes 
of the tuning and remains in the instant ground state if 
the system began from a ground state at t — > — oo; that 
is to say, the system undergoes an adiabatic evolution. 
Therefore, as Damski claimed, it is naturally suggested 
that the inverse of the instant energy gap acts as a quan- 
tum mechanical equivalent of the relaxation time scale 
r introduced above, which indicates how much time the 
system needs to adjust to new thermodynamic conditions 
|18j . As the direct results of the dynamics, the transition 
probabilities are suggested to identify the density of topo- 
logical defects - the main concern in the KZM. Then, as a 
direct suggestion, the other two fundamental parameters 
in the KZM, the relative temperature e and the quench 
time scale as tq, can be defined with At/ujQ and luq/A, 
respectively where ujq = 1/tq. 

In the KZM, all the predictions of the density of the 
topological defects can be made with the knowledge of 
the frozen-out time t, which satisfies Zurek's equation. 
According to the analogy presented above, this result 
means that the transition probability in the LZ model can 
be obtained from the knowledge of t. With simple mathe- 
matical calculations, Damski presented two cases to show 
that the transition probability, which corresponds to the 
density of the topological defects, can be predicted di- 
rectly and perfectly from the modified version of Zurek's 
equation by introducing a free parameter a to the equa- 
tion. Considering the case when the time evolution starts 
from a ground state at the anti-crossing center, the den- 
sity of the topological defects is suggested to be [18] , 

V n = 0.5(1 - ^l-2/V(x a )), (2) 

where x a — and V(x a ) = x 2 a + x a \/x% + 4 + 2. 

Due to the simple form of the Hamiltonian, the LZ 
model is experimentally obtainable and the transition has 
been observed in many quantum systems, from Rydeberg 
atoms [3011111] to quantum dots contacts [32], and recently 
extended to mesoscopic superconducting Josephson de- 
vices [33], ultracold molecules [33], optical lattices [35] 
and the NV center [3B] . There are no technological chal- 
lenges in controlling the level structure in LZ dynamics 
in most of these systems and more complicated dynam- 
ics based on the LZ model have been presented (see [37] 



for a review). Therefore, on the technological side, the 
verification of Eq[l] and the observation of the KZM in 
such a two-level dynamical system are possible. 

B. Characterization of the interferometer 

Two of the most important parameters of an interfer- 
ometer are the visibility and stability. Ultra-high visibil- 
ity interferometers have been used in various fundamen- 
tal physics, such as the laser interferometer gravitational- 
wave observatory [35]. Usually, the visibility is a direct 
indicator of the operation fidelity. The stability is an 
indicator of decoherence; a higher stability means less 
decoherence before the operations are finished. Here, we 
test our setup by measuring both the visibility and the 
stability, which are shown in Figure|4] The integration 
time of a single measurement is 1 s, and we monitor the 
phase stability for more than one hour. With the help of 
spectrum and spatial filter, the visibility is higher than 
0.998 for each stage. Then we connect the stages one 
by one with a lrh long single-mode fiber. The measured 
visibility decreases slowly to 0.975 ± 0.008 as the stages 
increase to nine, which fits the theory prediction of 0.998 9 
very well considering the error. We calculate the stan- 
dard deviation to investigate the stability in one hour. 
From the results shown in FigureQb), we can see that 
the stability worsens as the number of stages increases. 
However, due to the intrinsically stable feature of the 
Sagnac type architecture, even a nine-stage device is suf- 
ficiently stable to conduct the experiment. This finding 
is also confirmed with the x matrix ((c) and (d) in Fig|4|, 
which is reconstructed by the quantum process tomogra- 
phy. For an ideal case, the matrix equals I. Here, the 
measured fidelity is 0.994 ± 0.011. 

C. Observation of the evolution in LZ model 

With the help of this nine-stage interferometer, we per- 
form the simulation of the evolution of the LZ model, and 
the results are shown in Figure{5] Here, we choose the 
parameters of A = 5 and r = 2 for each step. Due to 
the limitation on the number of stages, one and a half 
periods in the evolution of the LZ model can be observed 
clearly. The points agree with the theoretical predictions 
very well. 
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